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Abstract 

After a brief historical survey, the paper introduces the notion of entropic 
model sets (cut and project sets), and, more generally, the notion of diffrac- 
tive point sets with entropy. Such sets may be thought of as generalizations of 
lattice gases. We show that taking the site occupation of a model set stochas- 
tically results, with probabilistic certainty, in well-defined diffractive properties 
augmented by a constant diffuse background. We discuss both the case of in- 
dependent, but identically distributed (i.i.d.) random variables and that of 
independent, but different (i.e., site dependent) random variables. Several ex- 
amples are shown. 



Introduction 

Diffraction is one of the most important ways of identifying long-range order in math- 
ematical and physical structures. In this paper, we look at the effects on diffraction 
that occur in certain periodic and quasiperiodic point sets when the occupation of the 
point sites is taken stochastically rather than deterministically, with independence 
between the different sites. Under fairly mild assumptions, which are certainly valid 
for lattices and model sets, we show that the effect is simply one of scaling down the 
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diffraction pattern by a constant factor and adding in a constant background. In the 
case of lattices, this type of phenomenon is well known ||, 39|. What is new here is 
that it remains true for a large class of non-periodic structures (Theorems 1 and 2) 
and, in particular, for all regular model sets (also called cut and project sets). 

The results are, on the one hand, a suitable reminder of the difficulty of interpreting 
the meaning of diffract ivity, and, on the other hand, of the robustness of diffractivity 
under certain deformations and modifications of the underlying set. 

It might be interesting to quickly review the history of aperiodic order and take a 
look at the reasons why stochastic forms of aperiodic structures seem to be a natural 
extension beyond the world of strict perfection. In the early eigthies, a new type 
of ordered state was found, both experimentally Q and theoretically fjl| [0| . 



These discoveries, made independently of one another, created an enormous amount 
of scientific activity because the new ordered states, quickly dubbed quasicrystals, 
had properties previously thought to be incompatible with one another; namely, long- 
range orientational order, strong enough to produce sharp diffraction images, and 
at the same time non-crystallographic symmetries such as icosahedral []30], [HJ or 
twelvefold symmetry |I4| . 

Since nothing has ever been discovered for the first time, one might expect precur- 
sors of this, and this is indeed correct. Clearly, Penrose's famous tiling of the plane 
with fivefold symmetry was important, particularly when coupled with de Bruijn's al- 
gebraic analysis that showed that it was also diffract iveQ. Ammann investigated this 
further and also found the matching icosahedral tiling made from two rhombohedra. 
However, this was not generally known in the physics community (a brief remark can 
be found in p^j), and has never been published by him. Even his planar results were 
published only much later ]T[|, though some results are contained in [10]. 

In fact, the history goes back quite a bit before this. In the late thirties of this 
century, Kowalewski |HJ investigated the possibilities of filling Kepler's triacontahe- 
dron with the two rhombohedra mentioned above, but apparently did not realize the 
possibility of filling the entire space with them. 

Kepler himself was very much interested in space fillings in his time, and his famous 
plate of planar tilings (resp. seeds of them) shows a considerable patch of a tiling with 
pentagons, decagons and rhombi, see the first plate in fTD |. In modern terminology QJ, 
this would be in the same MLD-class as the famous Penrose tiling, in the sense that 
there is a derivation rule with a radius much smaller than the patch shown such that, 
on the size displayed, there is no way to tell the two tilings apart. This is probably 



If not explicitly specified, the proper references are obvious by the names given, and can be 
found in M. 
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not an accident because Kepler was well aware of the problems of space fillings, as 
were other people before him, such as Diirer who constructed a mechanism to create 
a fivefold twin made from pentagons and rhombi. Generally, the investigation of 
geometric form was well on the way. Diirer's polyhedron in his "melancholia" has 
puzzled generations of scientists and art historians - with a really promising solution 
being found only very recently by Hans-Ude Nissen . 

Coming back to this century, the development of the theory of incommensurate 
structures by the Nijmegen group and the new developments in the theory of qua- 
sicrystals showed that there is a lot more to the geometry and symmetry of the solid 
state than anticipated by ordinary school knowledge. Probably the most puzzling 
aspect in the beginning was the combination of perfect diffractivity (in the sense of 
a Bragg spectrum) with non-crystallographic symmetry. But again, the final expla- 
nation, in terms of the projection method [I1J, had a precursor, this time in pure 
mathematics. 

Harald Bohr, the younger brother of Niels, developed the theory of quasi- and 
almost periodic functions in the twenties. The basic idea was to describe non- but 
quasi-periodic functions as sections through periodic functions in higher dimensions. 
In this sense, the cut-and-project method owes a lot to his ideas. In the late 1960's, 
Yves Meyer studied the harmonic analysis of point sets in the context of algebraic 
numbers. In the process he rediscovered cut and project sets (here called model 
sets), though now in the much wider setting of locally compact Abelian groups, and 



introduced a very important class of ordered point sets, now called Meyer sets p4 
(see later for one characterization of Meyer sets). 

After this historically motivated introduction, let us come to the aim of this article. 
Although all the above mentioned connections might indicate that the (quasi-) crys- 
talline world is perfect, in the sense that the alloys displaying such diffraction images 
are, reality tells us nowadays that this is not so |27|, |IE|| . In fact, quite early it was 
pointed out by Elser that, in order to explain the stability of such alloys, one might 
need an entropic side of the picture, an idea that led the Cornell group to develop the 
idea of a random tiling. From a more mathematical point of view, this is not fully 
satisfying because quite a number of questions concerning the diffractivity, and even 
the well-definedness of some of the ensembles, are still unanswered. 

This is the point we want to consider and start to develop. However, we will not 
adopt the random tiling picture here, because it seems not yet fully in reach for a rig- 
orous treatment (compare and references therein for some recent developments). 
Instead, as an intermediate step, we shall rather consider a setup of ideal model sets, 
or even more general diffractive point sets, that are coupled to stochastic processes to 
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thin them out. This way, we can introduce some randomness into the picture, even 
with positive entropy density, without losing control of the diffractivity. This can be 
seen as a generalization of the diffraction theory of lattice gases which belongs to the 
standard body of literature, see 0, |, [39|] and references therein. 

In doing so, we will actually arrive at results that exactly meet the expectation, 
but with the extra benefit of providing proofs for them, i.e. making a good deal of 
folklore rigorous in this way. It will turn out that there is a natural extension of the 
diffraction theory of lattices gases beyond lattices, provided one uses an approach that 
avoids techniques based upon translation invariance. Although the actual methods 
and results employed from probability theory and ergodic theory are pretty standard 
in mathematics, they are much less familiar to physicists. 

Let us now briefly sketch how the article is organized. We start with a section on the 
diffraction of lattice gases (without interaction) and its connection to entropy density. 
This gives us the opportunity to review some well-known results in a different setting 
that matches the generalizations derived later. We hope that the reader can adjust 
to our approach that way without too much pain. This is followed by our general 
setup, where we introduce a rather general type of point sets which are accessible 
with our methods. Common examples such as model sets (or cut and project sets) 
are contained as special cases. 

The remainder of the article is then devoted to the diffraction of point sets with 
independent stochastic occupation of sites. First, the focus is on the situation of 
independent, but identically distributed (i.i.d.) random variables, the case most fre- 
quently studied. Theorem 1 gives the result for this case. This is illustrated by some 
examples, and model sets in particular. 

More general, and less obvious, is the treatment of independent, but not necessarily 
identically distributed, random variables, which leads to Theorem 2. Among the 
applications are weighted model sets and their stochastic counterpart, and, more 
specifically, weighted model sets where the weights are determined by a so-called 
invariant density @, ||. This way, we are able to keep certain aspects of point and 
inflation symmetries. We believe that this application is of particular value in the 
discussion of perfect versus random tiling order, as it really is a first step of an 
intermediate picture. 

Our concluding remarks try to relate the results to other investigations and to 
point towards the next steps that should be taken. 
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Diffraction from a lattice gas 

In order to keep things simple, and to familiarize the reader with our approach, we 
start with the description of the lattice situation and give proper definitions for the 
general setup later. Let T be a lattice in W 1 , i.e. a discrete Abelian subgroup of R™ 
such that IR n /r is compact P3| |. Equivalently, there are n linearly independent vectors 
bi, . . . ,b n , called the basis vectors of T, so that T = 1b\ © ■ ■ ■ © Zb n . Since we will 
be talking about Fourier transforms, we will also need the dual (or reciprocal) lattice 
of r, defined by 

T° := {x G R n | x ■ y G Z for all y G T} (1) 

where x ■ y denotes the Euclidean scalar product. 

Next, define Dirac's comb as the characteristic distribution 



$> (2) 



uj = 0J T := 

on T, where 5 X is Dirac's distribution at point x, i.e. 

(6 x ,<f>) := (f>(x) (3) 
for all test functions <fi. In particular, one gets 

(u,,0) = (4) 

which is well defined for all rapidly decreasing functions (Schwartz functions), hence 
u r is a tempered distribution |28| . 

To deal with diffraction, we need the corresponding autocorrelation distribution, 
7 W , of T, also called its Patterson function (although it is a distribution^). With the 
abbreviation 

r r := rn B r (o) = {x g r | |x| < r} , (5) 

7 W can be defined and calculated as follows 

where d is the density of T, i.e. the number of lattice points per unit volume. 



2 It would actually be slightly more appropriate to adopt the setup of measure theory, where 7^ 
would represent a tempered measure, sec |l2j for this complementary approach. 
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By the Fourier transform of a Schwartz function we mean 

i(k) := I e- 27Tik - x <j)(x)dx (7) 



which is again a Schwartz function 28]. The inverse operation is given by 



^(x) = / e 2nlx - k ip(k)dk . (8) 

This definition results in the usual properties, such as = and ip = ip. The 
convolution theorem takes the form 0i * 2 = 0i • 02 where convolution is defined by 



<Pi*<p2{x) = / <f)i{x - y)<f) 2 {y)dy . (9) 

Finally, the matching definition of the Fourier transform of a distribution T is 

(T,0) := (T,}) (10) 

for all Schwartz functions 0, as usual. 

Now, the diffraction by the lattice T is described by the Fourier transform of its 
autocorrelation 0, and in this case we have 7 r = d ■ u) r . To calculate the latter, we 



employ Poisson's summation formula for tempered distributions, cf. p. 254 of [[33 
which reads 

Y,6 X = d-J2$k (11) 

xer fcer° 

and can easily be proved from the corresponding Poisson summation formula for 
Schwartz functions. So we get 

% = d 2 ■ h (12) 

fcgr° 

which is the well-known result that the diffraction from point scatterers of equal 
strength on the points of a lattice is a pure point measure, consisting of periodically 
placed point measures on the dual lattice. Note that the strength of the peak at k — 
is d 2 , the square of the density of T, as it must be. 

Let us now move on to the corresponding lattice gas, i.e. to the point set obtained 
from T by removing points from it stochastically. To describe this, we define a new 
measure u s (with suffix s for stochastic) by 



^8 = 



J2v(x)8 x (13) 
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where 77 (x) is a random variable at site x which takes only the values and 1, meaning 
empty or occupied. We assume that these random variables are independent of one 
another and identically distributed, i.e. they constitute a (countable) family of i.i.d. 
random variables || [36[. In fact, we parameterize the probability P that rj(x) takes 
the value 1 by a number < p < 1, i.e. 

P({ V (x) = l}) = p , P({ v ( x ) = 0}) = 1-p. (14) 

With this definition, u s of ([13|) describes scatterers on the sites of a lattice, each single 
site being occupied with the same independent probability p. Clearly, the mean value 
of each random variable is m 1 = p, the second moment is m 2 = p and the variance is 
thus v = m 2 — (m x ) 2 = p{l — p). 

Refering now to the strong law of large numbers || [36| , we can deduce (details will 
be given below in a more general context) that, almost surely, 

lim mu Yl V(x)v(x ~y) = d-p(p+(l- p)5 yfi ) , (15) 

i.e. the limit exists and equals the right hand side with probability one. Here, 5 a ^ 
denotes Kronecker's delta. With this expression, one can calculate the new autocor- 
relation to be, with probability one, 

lu s = P 2 -Icu + d- Pi 1 -p)-Sq. (16) 
Fourier's transformation then gives 

% s = P 2 -% + d-p(l-p) , (17) 

i.e. we retrieve the point part of the full lattice case, reduced by a factor of p 2 as it 
should according to the reduced density of points, plus a constant diffuse background 
which is the absolutely continuous part of the diffraction. This term is related to the 
entropy density s = s(p), 

s(p) = -plog(p) - (1-p) log(l - p) , (18) 

which is a measure of the complexity of the ensemble of point sets we are actually 
dealing with in this example. The presence of entropy is usually connected to con- 
tinuous components in the diffraction spectrum, and one can clearly see that s(p) 
vanishes if, and only if, p — or p — 1, i.e. iff the random variables are sharp. Also, 
s(p) is maximal at p = 1/2, which corresponds to the value of p where the diffrac- 
tion (|i~7|) shows the largest amount of white noise. Finally, in this case, we have an 
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essentially unique relationship between the entropy and the background intensity, up 
to the symmetry p <-> (1 — p). 

Note that the continuous part vanishes if and only if the variance of our random 
variable vanishes, i.e. if and only if p = or p = 1. Also, the pure point part vanishes 
if and only if the mean of the random variable vanishes, i.e. if and only if p = 0. We 
shall meet this situation again later in a much more general context. 

General setup 

The previous section should have served to get a feeling of what we are after, and 
to introduce the type of notation we wish to apply. Let us now develop the theory 
with more precision and in more generality. In what follows, we will only consider 
uniformly discrete point sets A C M n here, i.e. points sets with the property that 
there is a positive radius e such that each point x G A can be surrounded by an 
open ball of radius e that does not contain any point from A other than x. With 
this assumption, the corresponding Dirac comb uj a defines a translationally bounded 
distribution (i.e., for each compact set K, there is a constant ck so that for all x G M n , 
( u a^Xk+x) — c Ki where Xs denotes the characteristic function of a set S). This is 
sufficient, though certainly not necessary, to make uj a a tempered distribution. As 
before, we write A r for the intersection A D B r (0). Now, we have to tie this together 
with Fourier analysis. 

Definition 1 Let A be a uniformly discrete point set of (existing) natural density 
d > 0, and let u = u A = XLgA ^ ^ e ^ s Dirac comb. We say that A has a natural 
autocorrelation if 



exists as a limit in the weak topology (i.e., as a limit of tempered distributions), and 
thus is a tempered distribution. Then, 7^ is called an autocorrelation distribution or 
simply an autocorrelation of A. 

Note that, if this situation applies, then 7^ is translation bounded (7^ inherits this 
property from uj) and is a distribution of positive type, compare |12|]. Let us briefly 
comment on the more general setup in terms of measures. If we only had existence 
of an autocorrelation as a measure in the vague topology, translation boundedness 
would guarantee that it is actually tempered - so, in our context of uniformly discrete 
sets, the restriction to tempered distributions is reasonable. 




(19) 
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Note that existence of an autocorrelation, as we have defined it here, is specific to 
the type of region (in this case balls, as implied above by the attribute "natural" ) over 
which we compute our averages. Replacing balls by other (convex) objects, centered 
at 0, the limit r —>■ oo (with r the radius of inscribed balls, say) might give a different 
answer or might not even exist. For the purposes of this article we do not need to 
deal with limiting processes over more than one type of shape at the same time, so 
we will phrase the arguments in terms of convergence based on sequences of balls. 
The arguments for sequences based on other shapes work in the same way. However, 
for other purposes it is important to specify uniqueness of the autocorrelation. For 
example, this can be done as follows: 

Let (C n ) ng N be a family of convex bodies, centered at 0, with the properties that, 
as n — > oo, the radius of the maximal inscribed balls tends to oo and the quotient of 
the radii of circum- and inscribed balls is bounded. If 

7. := lim -J— Yl 5 *-v ( 2 °) 

exists for each such sequence as a tempered distribution and is unique, we say that A 
has a unique autocorrelation, 7 W . 

In any case, a point set with a natural autocorrelation has a positive measure 
7 W as its Fourier transform (due to the Bochner-Schwartz theorem |28fl). It is this 



measure that desribes the diffraction |9], [12[, and it is very natural that a positive 



measure shows up here: after all, diffraction is all about the amount of intensity 
scattered into a certain (measurable) reagion of space. Also, this positive measure 
can now, according to Lebesgue's decomposition theorem, uniquely be decomposed 
into an absolutely continuous part (ac), a singular continuous part (sc), and a pure 
point part (pp). The pure point part (usually called "Bragg part" in physics) will 
always contain a trivial term of the form d 2 ■ 5q where d is the (by assumption existing) 
density of A per unit volume. This motivates 

Definition 2 A point set A with autocorrelation r ) w is called diffractive (with respect 
to the convergence process adopted) if (l^pv ^ s non-trivial, i.e. contains Dirac distri- 
butions different from d 2 ■ 5 . A is called perfectly diffractive or pure point if has 
no continuous part at all. 

The simplest example of a perfectly diffractive point set is a lattice, where the 
statement follows from Poisson's summation formula. Another class of examples is 
given by regular model sets with sufficiently nice windows (see the last section for 
more on this) or by extensions of them to certain limit-periodic or limit-quasiperiodic 
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point sets usually described by means of inflation JJ7], 0. All these examples are not 
only uniformly discrete, but also relatively dense, so they are Delone sets. What is 
more, they are actually Meyer sets, i.e. not only are they Delone but they have the 
additional property that their difference set, A — A, is also Delone. Note, however, 
that the Delone property is not necessary for perfectly diffractive sets, as can be 
seen from the example of the set of visible points of a lattice || which has holes of 
arbitrary size (and this even with positive density) and is thus neither Delone nor a 
density deviation of one. Note that removing or adding points of density zero from 
a perfectly diffractive set does not change its autocorrelation, and the set thus stays 
perfectly diffractive. 

On the other hand, Meyer sets need not be perfectly diffractive, as can be seen from 
the union of 2Z with various subsets of 2Z + 1. This is always a Meyer set, but one 
can easily construct cases with continuous components (and positve entropy density). 
This indicates that the class of Meyer sets, or even Delone sets, and the class of 
perfectly diffractive sets are rather different, though they have some sets in common. 
In general, perfectly diffractive sets will not be Delone, and hence not Meyer. One 
interesting class of point sets in this context is that of uniformly discrete sets S with 
the extra property that S — S is Delone, or at least that S — S is closed and discrete. 
They are the ones we shall consider here. 



Point sets with independent stochastic occupation of sites 

In this Section, we will develop an appropriate generalization of the lattice gas (with 
i.i.d. random variables) to much more general point sets. From now on, let A be a 
uniformly discrete point set which has a natural autocorrelation. Let us also assume 
that A is of finite local complexity, i.e., that A := A — A is discrete and closed, compare 



21| for a detailed discussion in the context of Delone sets. Finite local complexity 
of a set A implies that, for every radius r > 0, there are, up to translations, only 
finitely many different configurations of points in a ball of radius r. In particular, A 
is uniformly discrete. This is so because A discrete forces G A to be isolated, so 
different points in A must have a uniform minimal distance from one another. 

If uj — J2xeA * s the Dirac comb of A, as usual, it is now certainly translation 
bounded, and the autocorrelation is given by Eq. ([19]). We then have 

7. = (21) 
zeA 
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where u(z) is the autocorrelation coefficient at z, defined by 

x,j/GA r 
x-y=z 

Here, A r := Afl-B r (0) as before, and V r := vol(_B r (0)). Note that, in order to establish 
the existence of an autocorrelation, it is sufficient to show the existence of the limits in 
(|22j), i.e. the existence of the coefficients, because A discrete then implies existence of 
the autocorrelation as a measure in the vague topology, and translation boundedness 
ensures temperedness, see |0| for further details. 



Let us now turn to a stochastic "lattice gas" version of A. It is defined by the 
characteristic distribution 

u s = ^2r](x)5 x , (23) 
where 77(2) is a family of i.i.d. random variables taking the values and 1, parame- 



terized as in Eq. ([14]), each with mean p and variance v = p(l — p). 

We first address the question of the existence of the corresponding stochastic au- 
tocorrelation. In analogy to Eq. (|22"|), we now have the coefficients 



u„{z) := lim -j- V r]{x)r](y) = lim -j- Y] r](x)r](x - z). (24) 

1 — >oo V r — r— >oo V r — 

x,yeA r x,x-zeA r 
x-y=z 

We will show that under mild assumptions these coefficients exist, at least in a 
probabilistic sense. In order to do this, we need to be able to decompose the sum 
involved in v s (z) into two parts, because the various terms in the sum of (|24"D are still 
random variables, but not necessarily independent ones any more. Fix z G M. n . Define 

S(z) := {x I x,x - z E A} (25) 

and its restricted version S(z,r), where the x, x — z appearing in the definition are 
required to lie in A r . We distribute the points of S(z) (and, by proper restriction, also 
those of S(z,r)) into two sets S(z)^ and S(zY 1 \ This may be done in an arbitrary 
fashion, subject only to the two conditions that 

(1) if x, x — z both lie in S(z), then they are not in the same S(zY' l \ and 

(2) the two sets S(z)^ have well-defined densities: 



u^(z) := lim — V 1. (26) 

r^oo V r 

xG5(z,r)( i ) 
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Evidently, ( fffi ) implies v(z) = is(°'(z) Let us say that the set A can be 

decoupled if, for every z G A, we can find such a partition. 

We note that there may be many ways of decoupling a set A. For lattices, for 
example, we can take each line of points x + Zz and distribute it into the two subsets 
according to whether the coefficient of z is even or odd. For aperiodic model sets (see 
definitions below), each of the sets {x \ x+z, x — kz ^ A ; x, x — z, . . . ,x—(k—l)z G A} 
is finite, of bounded length k and, for each k, the set of such strings has a definite 
density. We can place the points x, x — z, x — 2z, . . . alternately into S(z)^ and S(z)^ 
and again obtain sets with well-defined density this way; see |13j for a very similar 
approach to thermal fluctuations which establishes the usual form of the Debye- Waller 
factor for essentially the same kind of structures that we are dealing with here. 

The decoupling property is some kind of ergodicity assumption. It is certainly 
fulfilled for point sets with uniform frequencies of all finite patches (as is the case for 
usual model sets), but it is more general than this. In particular, it is still valid for 



objects such as the pinwheel tiling, compare the brief discussion in [Tj|. At present, 
we do not know any equivalent characterization simpler than that given above, which 
is very much designed for its (technical) purpose. 

Proposition 1 Let A be a point set of finite local complexity which has a natural 
autocorrelation. Assume further that the set A can be decoupled. Then each coefficient 
of the stochastic autocorrelation (i.e. the corresponding limit) exists with probability 
1, and is given by 

u s (z) = v{z) ■ (K) 2 + (m 2 - K) 2 )5 Z)0 ) (27) 

where m\[= p) is the common mean of the i.i.d. random variables r](x) and rri2 is 
their common second moment. 



PROOF: This is an application of the strong law of large numbers. Forming a sequence 
of random variables out of a family {i]{x)) x< za etc. is rather canonical. Since A is 
uniformly discrete, we number the rj(x) with x in finite sets A r for increasing r. Each 
such sequence, by the general assumptions made, is a sequence that conforms to the 
strong law of large numbers. 

Let us consider z = first. Here, the relevant random variable is actually r](x) 2 , 
with mean rri2, the second moment of rj(x). These variables are independent and, 
almost surely, 

lim — — y i](x) 2 = d ■ vrt2 

r^oo V r ^ — ' 
x—z£A 
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where d = z/(0) is the (existing) natural density of A. 

Next, let z 7^ 0, z G A, be arbitrary, but fixed. If v(z) = 0, also v s (z) = 0, and 
our assertion is trivial. So, assume v{z) > 0, which means that the density of points 
x G A, such that also x — z G A, exists and is positive. For each such x, r)(x)r)(x — z) 
is a random variable with mean (mi) 2 , where m\ = p is the (identical) mean of all 
random variables r)(y) involved. We now have to consider 



lim — i](x)i](x — 

r — *no V ' * 



r— >oo V r 

x£A r 
x—z£A r 

This sum has only non-negative terms and decomposes as two sums: 

v (z) = lim — i](x)i](x — z) + lim — i](x)i](x — z) . 

r^oo V r r^oo V r — ' 

i£S( Z ,r)(°) x£S(z,r)W 

Now each of the two sums is an averaged sum over a set of independent random 
variables. Hence, by the strong law, we get almost sure convergence to 

u s (z) = v {Q \z)m\ + v {1 \z)m\ = v(z)m\ (28) 

because the mean of each random variable r\{x)r\{x — z) is m 2 . Together with the first 
step, this establishes our claim. □ 

The autocorrelation 7^ of u s is defined as the distribution whose value on any test 
function <f) is 

(7uv0) = lim Tr E v(x)v(y)H x - y) ■ ( 29 ) 

x,yeA r 

Although we know this already from the above abstract arguments, it might be in- 
structive to check explicitly that 7^ is indeed a tempered distribution, at least in the 
sense of almost sure convergence. First, let be a C°° function of compact support, 
lying in the ball B s (0) of radius s. Then 



/— X) ' * \' 

zeA 



( \ 

. x,yeA r 

\x-y=z / 



(z) . (30) 



This limit exists because in reality the outer sum is over the finite set A s and for 
r » s we have A s C A r — A r . Thus, as r — > 00, the sum converges, almost surely, to 

(31) 
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Now if (j) 6 «S, the space of Schwartz functions, and {(pi} is a sequence of C°°- 
functions of compact support that converge to in the standard topology of S, then 

= (mO^i/^^^ + K-K) 2 )^)^^), (32) 

zgA zGA 

and the latter converges in i to 

(m x ) 2 ^ z/(z)0(z) + (m 2 - (m 1 ) 2 )z/(O)0(O), (33) 

zGA 

by our assumptions on the existence of the autocorrelation density of A. 
Let us summarize these findings as follows. 

Theorem 1 Let A be a point set of finite local complexity which has a natural auto- 
correlation and density d. Suppose that A can be decoupled. Then, the autocorrelation 
of A and that of its stochastic version are, with probability one, related by 

7o, s = (mi) 2 7 w + d(m 2 - (mi) 2 )5 - (34) 
As a consequence, their Fourier transforms fulfil 

% s = (mi) 2 i + ^K-(mi) 2 )- (35) 

So, the stochastic version has the same 'main' part of the diffraction, multiplied 
by a factor of (mi) 2 (hence vanishing if and only if the mean of the joint probabil- 
ity distribution is 0) plus an extra absolutely continuous part that is constant and 
represents the 'white noise' of the uncorrelated random processes. The constant is 
essentially given by the variance of the joint distribution, and thus this part vanishes 
if and only if the i.i.d. random variables are all sharp. The interpretation, and also 
the connection with the entropy density, is thus the same as in the lattice case, as 
expected. 

At this point, generalizations are rather obvious, and we just want to mention a 
few. First of all, it is by no means essential to restrict to the particular types of 
random variables that we have just discussed. Here we were motivated by the idea of 
a lattice gas and its generalization to uniformly discrete point sets, but we can also 
think of any other (non-negative) i.i.d. random variable with (existing) mean mi and 
second moment m 2 (so, the variance would be v = m 2 — (mi) 2 ). This does not change 
the result, and would correspond to a situation where we place, at each point x of the 
set A, a scatterer of random strength r](x). Again, we get the result of Theorem 1. 

Note also that at no point did we need to assume that 7^ was pure point. This is 
not necessary, indeed, and the result of Theorem 1 also applies to situations where 7^ 
is singular continuous, absolutely continuous, or of mixed type. This same situation 



is met in Hof's treatment of thermal fluctuations 13 
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Applications to lattices and model sets 

The obvious first application is to lattices. This results in a rigorous derivation of 
what we described in Section 2. The diffraction from a lattice gas, with i.i.d. random 
variables for the strength of the Dirac distributions at the lattice points, shows a point 
part that is the one from the lattice itself, reduced in intensity, plus a homogeneous 
diffuse background. 

Another application is to characteristic decorations on tilings^ that are obtained 
by a primitive substitution rule. Here, it was shown \22\ that the autocorrelation is 
unique, and convergence of its coefficients is even uniform. In general, the Fourier 
transform will not be pure point, see |37]] for a more detailed discussion. We also refer 
to jn| for a brief discussion of the decoupling property in situations without finite 
local complexity such as the pinwheel tilings of the plane. 

Lattice gas versions of model sets provide another class of examples, of rather 
recent interest. Recall that a model set [ESJ, is defined via projection onto R n 
of a lattice in some higher dimensional space, or, more generally, in some locally 
compact Abelian group. More precisely, it is assumed that G = R n x H is a locally 
compact Abelian group and that D is a lattice in G. Thus D is a discrete subgroup 
of G for which the quotient space G/D is compact. Further we assume that the 
projection vr 1 (D) of D into M n is injective and its projection tt 2 (D) into H is dense, 
where n 1 and 7r 2 denote the canonical projections. The resulting set is aperiodic, 
i.e. has no translational symmetries, if and only if ir 2 is injective on D. The most 
common examples take H = W" 1 for some m. In any case, define the composite map 
* '■= ^2 ° 71 i\d 1 '■ ^li-D) — ► H. Then for any set VL C H with nonempty interior and 
compact closure, we have the model set 

A = {x E n^D) | x* E CI}. (36) 

Provided that the boundary of Q has measure (with respect to the Haar measure 
/i of H), the density of such a set exists uniformly and is given by //(n)/vol(.D). Here, 
vol(-D) is the volume of any fundamental domain for D in G, the volume taken relative 
to the product measure on G derived from the Lebesgue measure on MJ 1 and the Haar 
measure /i on H [RT], |3"2"fl . Such a model set is a Meyer set, i.e. both A and A = A — A 
are Delone. Also, A is perfectly diffractive | 3~3]| , and the obvious lattice gas version of 
it, with i.i.d. random variables attached to each position, falls under our Theorem 1. 

A large number of well-known point sets can be interpreted in this setting, including 
the Fibonacci and many other chains, the vertex sets of various planar tilings (such as 

3 We call a point set of finite local complexity a characteristic decoration of a locally finite tiling 
if they are locally equivalent, i.e. if both objects represent the same MLD-class, see (fij for details. 
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the Ammann-Beenker, the Penrose, the Tubingen triangle tiling etc.) or of tilings in 
higher dimensions (such as the various icosahedral examples in 3D or the Elser-Sloane 
quasicrystal in 4D). But even decorations of the chair tiling and other limit-periodic 
and limit- quasiperio die structures fall under this class, see |§ for details. So, for all 
these cases, we have 

Corollary 1 If A is a model set as described above, it fulfils the conditions of Theorem 
1, and the diffraction of the stochastic versus the deterministic Dirac comb is, almost 
surely, given by Eq. ffidj). 

Beyond identical distribution 

So far, we have restricted our attention to the case of i.i.d. variables. We will now 
broaden our point of view to the situation where the random variables are still inde- 
pendent, but not necessarily identically distributed any more. Before we give proper 
definitions, let us have another look at the lattice gas. Above, we compared the deter- 
ministic Dirac comb uj = ^2 xer S x with the stochastic one, uj s = 'Yl, x( z T 'u{x)5 X) where 
rj(x) were i.i.d. random variables of common mean m l . This led to Theorem 1. 
Alternatively, consider now the deterministic, but weighted distribution 

u a ■= ^2™> X 8 X . (37) 

Clearly, uj a has an autocorrelation if uj itself does, and we get 

lu> a = K) 2 "7. (38) 
and the result of Theorem 1 may be restated as 

% s = % a + d-(m 2 -( mi ) 2 ) (39) 

which holds almost surely. 

This indicates how we have to generalize our previous findings properly. Let A 
again be a uniformly discrete set of finite local complexity (i.e. A = A — A discrete 
and closed), and suppose that A can be decoupled in the sense described above. 
Furthermore, let (r](x)) x& \ be a family of independent random variables with non- 
negative means m x (x) which are bounded from above and with bounded variances 
v(x), v(x) < c, say. Under these assumptions, this family conforms to the strong law 
of large numbers. This can be seen as follows. Let (r] m ) m( z?>s be any sequence made 
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from the random variables r](x), e.g. by numbering the points of A in balls of growing 
radius. We then obtain 

E^^E^— C(2) = =-<». ,40) 

m=l m=l 

The assertion now follows from Kolmogorov's criterion, see or Theorem 12.3], 
and we can continue to develop the appropriate analogue of Theorem 1. 
To this end, let us now compare the two distributions 

u a = ^m 1 (x)5 x , (41) 

xeA 

which may be thought of as a toy model for an arrangement of different atoms (hence 
the suffix a), and 

the former being deterministic and the latter probabilistic. We can now formulate the 
appropriate theorem for this situation. 

Theorem 2 Let A be a set of finite local complexity that can be decoupled, and let 
(v( x ))xeA be a family of independent random variables with non-negative means m x {x) 
(bounded from above) and bounded variances v(x) whose average is assumed to exist, 

v = lim — — - v(x) . (43) 

r-+oo AJ 

1 1 x£A r 

If oo a has a natural autocorrelation, in the sense we used this term above, then uo s also 
possesses, almost surely, a natural autocorrelation, namely, with d = dens(A), 

lu a = + (44) 

and its Fourier transform reads 

% g = % a + dv. (45) 

The proof is very similar to the one given above and need not be repeated. 

This theorem is a little less explicit than the previous one, and one can see the 
potential extra complication from the following simple example. Consider A = Z and 
independent random variables 77(771) with values in {0, 1} and parametrization 

r p, m even, 

P({„(m) = l}) = {I mQdd ; (46) 
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where < p, q < 1. Here, u a = q^2 xeZ S x + (p — q)^2 xe2 z^x which clearly has a 
unique autocorrelation 7^ , with Fourier transform 

(47) 

So, the diffraction spectrum depends on the values of p and q, and the second term 
on the right hand side vanishes for p = q. 

The corresponding stochastic version, cu s , reflects this and produces the same point 
diffraction, plus a constant diffuse background ("white noise"), i.e. by application of 
Theorem 2 we have, almost surely, 

% s = i Ua + \w-p) + Q(i-q)]- (48) 

The entropy density of this little example is immediate: 

s = -^[P i og(p) + (l-p)log{l-p) + qlog{q) + (l-q)\og(l-q)]. (49) 

Note that, in general, a perfectly diffractive point set A together with a family 
of independent random variables r](x) is not enough to apply Theorem 2; we really 
have to know that not only 7^ but also 7^ exists. This is a rather subtle (and non- 
constructive) set of conditions upon the means of the random variables. There is one 
situation where we can escape this extra complication: if the random variables are 
distributed statistically, i.e. in such a way that their means m 1 (x) are themselves the 
result of a stationary Bernoulli process, we are back to the situation of Theorem 1, 
which may then be applied with 

m i = lim TT1 m i( x ) ' ( 50 ) 



x£K r 



provided this limit exists. 



Further examples: weighted model sets 

Let us now come back to the situation of a model set A = A(fi), as described above. 
Assume that we have a family of independent random variables parametrized by the 
points x* of the window Q. Suppose that p(x*) is a continuous function on Q, with 
values in [0,1]. Then, uj a = J2 x &aP( x *)^ is perfectly diffractive (this follows from 
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a slight modification of the arguments given in |T2| by means of an application of 
Weierstrass' approximation theorem). Explicitly, one has 



where D° is the dual lattice of D and the amplitudes are given by |T2| 



where d denotes, as before, the density of the model set A. 
Let us turn to the stochastic counterpart 



U s = 



^2rj(x)5 x , (53) 



where rj(x) is the random variable that decides whether x is occupied or not. Let us 
define it as follows 

P({r ] ( x ) = l}) = p(x*) , P({r](x) = 0}) = I - p(x*) . (54) 

Observe that r](x) has mean p(x*) and variance p(x*)(l — p(x*)), the latter being 
bounded by 1/4. So, by Kolmogorov's criterion, this family of random variables 
conforms to the strong law of large numbers. 

Let us see whether the mean of the variances exists. We note first that we have 
the mean occupancy per point of A as 

V = lim 7~7 r Y] p(x*) . (55) 

r-»oo A r z — ' 
1 1 x£A r 

Due to the fact that the projection into internal space H is uniform and the fact that 



p is continuous, it is possible to use Weyl's theory of uniformly distributed sets |20 
to show that this limit indeed exists and is given by 

P = TTTy, I P(v)dn(y). (56) 

vol(il) J n 

In the same way, we can also calculate the averaged variance as 

1 



vol(fi) 



p{y) 2 d^i{y)-p 2 . (57) 

n 
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\p(y) \og(p(y)) + (1 - p(y)) log(l - p{y))\ d^{y) . (61) 



So, we can apply Theorem 2 and obtain, with probability one, 

% s = % a +dv. (58) 

The resulting set also has a positive entropy density Clearly, for a single point, 
this is 

s(x) = — p(x*) \og(p(x*)) — (1 — p(x*)) log(l — p(x*)) (59) 
and we would be interested in the quantity 

1 := lim TT1 Yl s ^ ' ( 60 ) 

r— >oo \A r \ * — ' 
1 1 xeAr 

provided this limit exists. Again, this follows from the uniform distribution of the 
points of a model set and, using Weyl's lemma, we obtain 

1 

S ~ "vol (ft) 

As a relevant example, let us consider a special function p(x*), namely one that 
reflects the inflation structure of a given model set and is related to the recently 
investigated invariant densities on them || H|. Assume that G = MJ 1 x M. m and 
suppose that A = {x G 7Ti(D) \ x* E Q}. We are interested here in the situation in 
which A admits self-similarities of the form 

t Q)V : x h-> Qx + v , (62) 

where Q is an inflational linear map, i.e. a rotation followed by a scalar inflation. We 
call such self-similarities Q-inflations. Remarkably, for fixed Q, the set 

T := {v G R n | t Q>v A C A} 

is itself a model set. In this situation, there is a unique absolutely continuous proba- 
bility measure p = Pq supported on Q which is invariant under the set of Q-inflations 
in the sense that 



P 



lim-L^.p, (63) 



v£T s 

where t* is the induced mapping in internal space and (t* ■ p)(y) '■= p((t*) _1 y). 

The corresponding stochastic model set with site occupancy probability P({r)(x) = 
1}), is likewise invariant in the sense that 

P({ V (x) = 1}) = lim £ P({ V (y) = 1}) . (64) 
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Such invariant densities in internal space are supported on the window and typically 
display bell-shaped form. We refer the reader to [§, §, |[] for more details on this and 
for various examples. 

Since this is a special case of the general situation met above, the measure 

u s ^2v(x)S x (65) 

is (almost surely) diffractive with the same pure point part u a , since 

lim — r]{x)r]{x — z) = lim — p(x*)p(x* — z*) . (66) 

i — >oo v ^ — ' r^oo ~\/ r ^ — ' 

It is an interesting feature of this situation that the probability distributions pqk , 
as k — » oo, tend towards the constant distribution on Q. 



Concluding remarks 

The analysis of diffraction from point sets with stochastic occupation of sites, or 
with random scattering strength on the sites, can be developed in a rather general 
setting which goes considerably beyond the lattice situation. It was the aim of this 
contribution to outline some of the methods needed. For related aspects, we also 
recommend Hof 's treatment of thermal fluctuations [|13[ . 



One concrete reason to look into this type of problem stems from the discussion of 



quasicrystalline order and the evidence of stochastic elements in it [15, 16]. Based upon 
the random tiling scenario, one would expect the *-image (lift) of a "real world" point 
set (e.g. one obtained from a tiling overlay of a high resolution electron micrograph) 
to show a Gaussian shaped distribution or at least a bell shaped curve with maybe 
a somewhat flatter centre - in contrast to the uniform distribution obtained from a 
perfect model set. 

Since such bell shaped distributions have been observed and appear to be rather 
typical |16| , it is an important question to what extent they really support the random 



tiling picture. In other words: are there alternatives to explain such profiles? One is 
provided by the stochastic occupation of a model set, if we start from an invariant 
density on the window that resembles such a bell curve, see [Q], and [|J for the specific 
example of the Penrose tiling and invariant densities attached to it. A *-image of 
a finite patch its stochastic point set realization would reproduce the bell shaped 
invariant density profiles. 

We do not claim that this is enough to establish this simple stochastic approach as 
a real alternative - there are various other objectives to be met, such as width of the 
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profile as a function of the patch size, structured diffuse scattering background (other 
than white noise), or maximization of entropy as a function of suitable parameters for 
the non-crystallographic phase. We believe, nevertheless, that there are interesting 
possibilities along the lines presented here, and it would be nice to find (solvable) 
examples with extra correlations that appear more realistic in the sense mentioned. 

One first step in this direction is the calculation of the diffraction of a stochastic 
version of Z with random variables that stem from a stationary ergodic Markov sys- 



tem, as is well known in the literature [39]. This results in a pure point part which 



is that of Z, reduced in intensity, plus an absolutely continuous background that now 
shows a structure, i.e. that is no longer white noise. We hope to report on proper 
generalizations of this scenario soon, and also on an extension to the diffraction theory 
of random tilings. 

Acknowledgements 

It is our pleasure to thank Martin Schlottmann for several clarifying discussions and 
Vladimir Rittenberg for helpful advise on the manuscript. This work was supported 
by the German Science Foundation (DFG) and the Natural Science and Engineering 
Research Council of Canada (NSERC). 



References 

[1] R. Ammann, B. Griinbaum and G. C. Shephard, "Aperiodic tiles", Discrete Comput. 
Geom. 8 (1992) 1-25. 

[2] M. Baake and R. V. Moody, "Self-similarities and invariant densities for model sets", 
in: Algebraic Methods and Theoretical Physics, ed. Y. St. Aubin, CRM Publications, 
Springer, New York (1999), in press. 

[3] M. Baake and R. V. Moody, "Invariant densities for quasicrystals" , preprint (1998), to 
appear in Directions in Mathematical Quasicrystals, eds. M. Baake and R. V. Moody, 
CRM monograph series, AMS, Rhode Island (1999), in preparation. 

[4] M. Baake and R. V. Moody, "Multi-component model sets and invariant densities", 
in: Aperiodic '97, eds. M. de Boissieu, J. L. Verger-Gaugry and R. Currat, World 
Scientific, Singapore (1998), in press. 



22 



[5] M. Baake, R. V. Moody and P. A. B. Pleasants, "Diffraction from visible lattice points 
and kth power free integers", preprint (1998), to appear in Discrete Comput. Geom. 

[6] M. Baake, R. V. Moody and M. Schlottmann, "Limit- (quasi-)periodic point sets as 
quasicrystals with p-adic internal spaces" J. Phys. A31 (1998) 5755-65. 

[7] M. Baake and M. Schlottmann, "Geometric aspects of tilings and equivalence con- 
cepts", in: Proceedings of the 5th International Conference on Quasicrystals, eds. C. 
Janot and R. Mosseri, World Scientific, Singapore (1995), pp. 15-21. 

[8] H. Bauer, Probability Theory, de Gruyter, Berlin (1996). 

[9] J. M. Cowley, Diffraction Physics, 3rd ed., North-Holland, Amsterdam (1995). 

[10] B. Griinbaum and G. C. Shephard, Tilings and Patterns, Freeman, New York (1987). 

[11] A. Guinier, X-ray Diffraction in Crystals, Imperfect Crystals and Amorphous Bodies, 
Freeman, San Francisco (1963); reprint: Dover, New York (1994). 

[12] A. Hof, "On diffraction by aperiodic structures", Commun. Math. Phys. 169 (1995) 
25-43; 

A. Hof, "Diffraction by aperiodic structures", in: The Mathematics of Long-Range 
Aperiodic Order, ed. R. V. Moody, NATO ASI Series C 489, Kluwer, Dordrecht (1995), 
pp. 239-68. 

[13] A. Hof, "Diffraction by aperiodic structures at high temperatures", J. Phys. A28 
(1995) 57-62. 

[14] T. Ishimasa, H.-U. Nissen and Y. Fukano, "New ordered state between crystalline and 
amorphous in Ni-Cr particles", Phys. Rev. Lett. 55 (1985) 511-3. 

[15] D. Joseph and M. Baake, "Boundary conditions, entropy and the signature of random 
tilings", J. Phys. A29 (1996) 6709-16. 

[16] D. Joseph, S. Ritsch and C. Beeli, "Distinguishing quasiperiodic from random order in 
high-resolution TEM images", Phys. Rev. B55 (1997) 8175-83. 

[17] G. Kowalewski, Der Keplersche Korper und andere Bauspiele, K. F. Koehlers Anti- 
quarium, Leipzig (1938). 

[18] P. Kramer, "Non-periodic central space fillings with icosahedral symmetry using copies 
of seven elementary cells", Acta cryst. A38 (1982) 257-64, reprinted in [38]. 



23 



P. Kramer and R. Neri, "On periodic and non-periodic space fillings of E m obtained by 
projection", Acta cryst. A40 (1984) 580-7; and: Acta cryst. A41 (1985) 619 (Erratum), 



reprinted in [38|. 



L. Kuipers and H. Niederreiter, Uniform Distribution of Sequences, Wiley, New York 
(1974). 

J. C. Lagarias, "Geometric models for quasicrystals. I. Delone sets of finite type", 
preprint (1997), to appear in Discrete Comput. Geom. 

W. F. Lunnon and P. A. B. Pleasants, "Quasicrystallographic tilings", J. Math, pures 
et appl. 66 (1987) 217-63. 

A. L. Mackay, "De nive quintangula: on the pentagonal showflake", Sou. Phys. Crys- 



tallogr. 26 (1981) 517-22, reprinted in ||. 



Y. Meyer, Algebraic Numbers and Harmonic Analysis, North-Holland, Amsterdam 
(1972). 

R. V. Moody, "Meyer sets and their duals" , in: The Mathematics of Long-Range Ape- 
riodic Order, ed. R. V. Moody, NATO ASI Series C 489, Kluwer, Dordrecht (1995), 
pp. 403-41. 

H.-U. Nissen, "The polyhedron in Diirer's melancholia" , in preparation. 

H.-U. Nissen and C. Beeli, "Electron microscopy of quasicrystals and the validity of 
the tiling approach", Int. J. Mod. Phys. B7 (1993) 1387-413. 

M. Reed and B. Simon, Methods of Modern Mathematical Physics I: Functional Anal- 
ysis, 2nd ed., Academic Press, San Diego (1980). 

C. Richard, M. Hoffe, J. Hermisson and M. Baake, "Random tilings - concepts and 
examples", J. Phys. A31 (1998) 6385-408. 

D. Shechtman, I. Blech, D. Gratias and J. Cahn, "Metallic phase with long-range 
orientational order and no translational symmetry", Phys. Rev. Lett. 53 (1984) 1951- 
3, reprinted in [38]. 



M. Schlottmann, Geometrische Eigenschaften quasiperiodischer Strukturen, Disserta- 
tion, Univ. Tubingen (1993). 

M. Schlottmann, "Cut-and-project sets in locally compact Abelian groups", in Qua- 
sicrystals and Discrete Geometry, ed. J. Patera, Field Institute Monographs, vol. 10, 
AMS, Rhode Island (1998), pp. 247-64. 



24 



[33] M. Schlottmann, "Generalized model sets and dynamical systems" , to appear in: Direc- 
tions in Mathematical Quasicrystals, eds. M. Baake and R. V. Moody, CRM monograph 
series, AMS, Rhode Island (1999), in preparation. 

[34] L. Schwartz, Theorie des Distributions, 2nd ed., Hermann, Paris (1966). 

[35] R. L. E. Schwarzenberger, N -dimensional Crystallography, Pitman, San Francisco 
(1980). 

[36] Ya. G. Sinai, Probability Theory, Springer, Berlin (1992). 

[37] B. Solomyak, "Dynamics of self-similar tilings", Ergod. Th. &; Dynam. Syst. 17 (1997) 
695-738. 

[38] P. J. Steinhardt and S. Ostlund (eds.), The Physics of Quasicrystals, World Scientific, 
Singapore (1987). 

[39] T. R. Welberry, "Diffuse X-ray scattering and models of disorder" , Rep. Prog. Phys. 
48 (1985) 1543-93. 



25 



